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Abstract. We extend the theory for laser cooling in a near-resonant optical lattice to include multiple 
excited hyperfine states. Simulations are performed treating the external degrees of freedom of the atom, 
i.e., position and momentum, classically, while the internal atomic states are treated quantum mechani- 
cally, allowing for arbitrary superpositions. Whereas theoretical treatments including only a single excited 
hyperfine state predict that the temperature should be a function of lattice depth only, except close to 
resonance, experiments have shown that the minimum temperature achieved depends also on the detuning 
from resonance of the lattice light. Our results resolve this discrepancy. 

PACS. 32.80.Pj Optical cooling of atoms; trapping - 03.65.Sq Semiclassical theories and applications 



1 Introduction 

Laser cooling is a generic name for a number of techniques 
that use laser light to cool atoms down to millikclvin or 
even microkelvin temperatures [1] . Today, laser cooling is 
used in numerous applications, for instance as one of the 
steps used in the process to create a Bosc-Einstcin con- 
densate [2], in atomic clocks [3], and other high- precision 
experiments using atoms. 

The most commonly used technique, Dopplcr cooling, 
has as its lower limit the Doppler temperature, which for 
most atoms is of the order 0.1 mK [T] . However, lower tem- 
peratures can be reached in near-resonant optical lattices, 
i.e., standing waves of laser light with periodical, spatially 
alternating polarisations [4j[5j - In these systems temper- 
atures approaching the recoil limit of a few fxK can be 
achieved [BJ. This result first came as a surprise, but was 
soon given a theoretical model in form of the so-called 
Sisyphus mechanism 0IE]. In this model a combination 
of spatially dependent optical pumping rates between the 
magnetic sublevels of the atom, together with the periodic 
potentials of the optical lattice, give rise to an effective 
friction which causes cooling of the atoms. 

Whereas the Sisyphus model successfully describes many 
of the qualitative features of laser cooling in optical lat- 
tices, it is not sufficient for a quantitative analysis. For 
this purpose a number of numerical techniques have been 
developed, e.g., semiclassical methods based on Fokker- 
Planck-like equations [9lfT0] , a band-structure model [11] , 
and quantum Monte Carlo simulations [T2] ■ (For a review 
see reference [5].) These theoretical techniques have gone a 
long way in reproducing experimental findings. However, 
unexplained features still remain. For instance, theoret- 
ical simulations have consistently given kinetic temper- 
atures of the atoms which are independent of the laser 



detuning from the atomic resonance, except very close to 
this resonance |13| . Experimental results show that this is 
largely true for large potential depths (large laser irradi- 
ances), where the kinetic temperature is a linear function 
of the potential depth. However, as the laser irradiance is 
lowered a minimum temperature is achieved, before the 
temperature starts to rapidly increase for even lower ir- 
radiances. The point of this minimum is often referred to 
as decrochage. Experiments have shown that, in contrast 
to theoretical predictions, the point of decrochage docs 
depend on detuning [UHfToTfTB] . 



Hitherto all theoretical simulations have used a sim- 
plified level structure of the atom. It has been assumed 
that the cooling process only depends on optical pumping 
via a single excited state. However, the excited state is 
really a manifold of several closely-lying hyperfine states. 
In caesium the excited state used has total angular mo- 
mentum F c = 5. However, separation between the F B = 5 
and F c = 4 excited states is only 48. IP (r being the nat- 
ural lincwidth), which is comparable to typical detunings 
in experiments. A recent experiment showed that even for 
detunings very close to the F c = 4 state, the proximity to 
this state did not seem to have any effect on the temper- 
ature |15j . This seems to underpin the assumption that 
this state can be neglected in simulations. Nevertheless, 
it is still possible that the F c = 4 state is important for 
the cooling process close to decrochage. In this paper this 
possibility is investigated using semiclassical simulations 
including both the F c = 5 and F c = 4 states of the hyper- 
fine manifold. 
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2 Method 

In an earlier publication we developed a novel scmiclassi- 
cal method for Sisyphus cooling [10] , and showed that this 
method gives excellent agreement with the fully quantum- 
mechanical method |12j . In the semiclassical method the 
external degrees of freedom, i.e., position and momentum, 
are treated as simultaneously well-defined classical vari- 
ables. The internal degree of freedom, i.e., the magnetic 
substate, is on the other hand treated fully quantum me- 
chanically, allowing for arbitrary superpositions. In this 
way we are able to generalize the simplified F g = 1/2 — > 
F c = 3/2 model to realistic angular momenta, F g = 4 — > 
F c = 5 for Cs, while retaining an excellent agreement with 
fully quantum-mechanical simulations. This is in contrast 
to, e.g., the treatment in reference [9] where the internal 
states were projected onto an adiabatic basis. Our method 
automatically includes all couplings between adiabatic (or 
diabatic) states, which were neglected in reference [9]. In- 
clusion of these couplings has been showed to be crucial 
for good agreement with fully quantum-mechanical simu- 
lations [TO] . 

Having established the validity of our method, we now 
continue onto more detailed investigations of the cooling 
process. We first extend our method to include two excited 
hyperfine states. The generalized optical Bloch equations 
for an atom of mass m with a ground state g and two 
excited states el and e2 are 



da 
~dt 



P 

2m 



H A + Val,<t 



dcr 



(1) 



sp 



Here a is the density matrix of the atom, including the 
ground state g and both excited states el and e2, 



Sg "'gel a gc2 
Cclg Cold (7 c io2 
<7o2g Ce2el 0"e2e2 , 



(2) 



elements G e are products of the appropriate transition 
dipoles d and the positive frequency component E + of 
the laser field 



HG (r) =d-E+(r), 
E+(r) = EoS(t), 



(5) 
(6) 



with £ the (position-dependent) polarisation vector. In the 
basis of the magnetic substates, M e = —F e , —F e +1, . . . ,F C 
(c = el or e2), the transition dipole is the product of a 
reduced matrix clement and a Clebsch-Gordan coefficient. 



d c 
d c 

dl 



m\d\ 



f d° e e - 



d: 



(F g lM g q\F e M e ), 



(7) 
(8) 
(9) 



and e q are the usual spherical polarisation vectors. The po- 
larisation vector ^ is chosen as the one-dimensional lin_Llin 
laser configuration [5], 



£(z) = cos(fcz)e_i — isin(fcz)e + i, 



(10) 



with k the wave vector of the laser. The final term in equa- 
tion ([1]) describes the transfer of populations and damping 
of coherences due to spontaneous emission 



da 
dt 



7g 



Celg-fel/2 — Cclcl^cl 



"Cgel-fel/2 — <7 g e2-f n e2/2 



sp \ -cr o2g r e 2/2 









— Ce2c2^c2 



(11) 



where r c is the partial width of the excited state e for 
decay to the ground state g. (Where decay to hyperfine 
states other than g are possible, most experimental set- 
ups include a repumper laser which brings the atom back 
to the excited state.) Here 7 gg includes both the recoil 
and the probabilities of populating different ground states 
after a cycle of optical pumping through cither of the two 
excited states, 



where each cry is a submatrix with rows and columns cor- 
responding to the different magnetic sublevcls. The Hamil- 
tonian part of the evolution is determined by the kinetic 
term, the atomic internal Hamiltonian Ha and the atom- 
laser interaction Val- Setting the zero of energy at the 
ground state, the atomic Hamiltonian is just 



'0 
H A = | Hlo Fl!1 

,0 fiWF e2 



(3) 



with huF B the energy of the excited states. The atom- laser 
interaction takes the form 



o gU^gUvY 

Gci(r) 
t G e2 (r) 



(4) 



where G e represents the simultaneous absorption of a pho- 
ton and excitation of the atom in the excited level e, while 
G\ represents the inverse emission process. The matrix 



7gg = ^ / £[d ol • e]t e — r a clcl e lKr d cl • e 
+ '^ljdn K J2 [dc2 • e]t e - iK r a c2c2 e lKr d c2 • e. (12) 

In most cases of interest, the irradiancc of the lasers is 
sufficiently low that the population of the excited states is 
much smaller than the ground-state population, i.e., the 
transitions are far from saturation. This condition can be 
expressed in terms of the saturation parameters s e i and 
s e 2 as 

Ql/2 



Al + ri/4 



< 1, 



(13) 



where fi c = — (F c \\d\\F g ) E / h is the Rabi frequency of the 
transition and A e the detuning from the excited state. In 
the low saturation limit the excited states will rapidly ad- 
just to any change in the ground state density matrix. 
For time scales relevant to the evolution of the ground 
state, the excited state density matrices can be expressed 
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/gg. Through the process of adiabatic elim- 



as functions of <r g 
ination of the excited states we then arrive at an effective 
equation for <7 gg only, 



The first term above is the force arising from the second- 
order light-shift potential, while the second term is the 
radiation pressure. The diffusion coefficient is given by 



dOgg = J_ 

dt m 
3H 



r' 



Hcff,cr g g -^2^-{A c (z),a gg } 



r'h 2 h 2 



„ , , T OqQ «- 



9=0, ±1 



ES/^E^^^W' ( M ) + E^ E Tr{B° q '(z)B?(zMt)}, (23) 



where 



bi(z)= dt-r( 



dp • 



(15) 



Here the effective Hamiltonian describing the conservative 
part of the evolution is given by 

P 2 



H„ 



2m 



The matrix 



A c (z)= Idle (z)] \d e -£(z) 



(16) 



(17) 



describes the coherent couplings and light shifts of the 
magnetic substatcs. We have also used the conventional 
notation 

k=^, n = ^. as) 

In our simulations a scmiclassical approximation to 
equation (|14[) is used. This approximation is derived by 
first rewriting equation (|14p in terms of the Wigner dis- 
tribution, which is then Taylor-expanded to second or- 
der in p. For more details see [10]. This results in a set 
of coupled Fokker-Planck-like equations for the popula- 
tions of and coherences between the magnetic substates. 
Hence, while position and momentum are treated as clas- 
sical variables, the internal states of the atom are treated 
fully quantum mechanically. Finally, the evolution equa- 
tions are converted into a Langevin form. That is, each 
atom is assigned a time-dependent position z(t) and mo- 
mentum p(t) (see, e.g., [IB]), which follow the classical 
equations 



P 

z = — , 
m 



p = f(t)+ V (t). 



(19) 
(20) 



Here, f(t) is a conservative force and r){t) is a diffusive 
force with the properties 

<»?(*)> = 0, ( V (t)r 1 (t')) = 2D(t)S(t-t'), (21) 

where (■) stands for a time average and D is a diffusion 
coefficient. The forces are calculated as a trace over the 
magnetic sublevcls, where the internal state of an atom is 
represented by a density matrix w(t). The force is given 

by 



iEf E Tr{[B e q (z)Bf(z) 
e 9=0, ±1 

-B q '(z)B?(z)]w(t)}. 



(22) 



9=0, ±1 

with &;r the wave vector of the emitted photon (we ne- 
glect the difference in energy of the photons emitted from 
the two excited states). The first term arises from the re- 
coil from photons spontaneously emitted in random direc- 
tions, while the second term is connected to fluctuations 
in the radiation pressure. The evolution equation for the 
internal-state density matrix is 

w(t) =E {iA' c [w(t),A c (z)} - ^{ w (t),A e (z)} 



II E B?(z)w(t)B q (t 

9=0, ±1 



(24) 



Finally, we consider the specific case of the F g = 4 and 
F i = 4, F C 2 = 5 states of caesium. Using that (4||d||4) = 
v / 77l2(5||d||4) [IH] we find that H 4 = (7/12) Q 5 . We also 
have I5 = r, where r is the natural lincwidth, while for 
the the excited state F e = 4 the partial width for decay to 
the F g — 4 ground state is (7/12) r. The energy separation 
between the excited states gives A4 = 48 Ar + A5. From 
these relations the respective saturation parameters (|13|) 
and r^, A' e , equations p8]) . can be derived. 



3 Results 

We have performed simulations using 5000 atoms for de- 
tunings A 5 = -10r, -20r, -30r, -40r and a range 
of different potential depths (directly proportional to the 
parameter A' 5 ). In all simulations the initial temperature 
was 10 /LtK and the time step was dt = 0.001//^. The 
simulations were iterated until the second moment of the 
momentum distribution had stabilized. It should, however, 
be noted that since simulations are necessarily performed 
using a finite number of atoms the moments of the mo- 
mentum distribution will still fluctuate over time irrespec- 
tively of the number of iterations. The typical size of these 
fluctuations were A(p 2 ) < p^ (where is the recoil mo- 
mentum), but grow for potential depths below decrochage. 
Depending on detuning and potential depth the number of 
iterations required for convergence varied between 100000 
and 600000. The stability of the results were checked with 
larger numbers of iterations and with shorter step sizes. It 
was found that the step size 0.01//^, which was used in 
reference [10], while working well at A$ = — 10T was too 
coarse for larger dctunings. 

Our main results are displayed in figure [TJ where the 
one-dimensional temperature T, defined as k^T '/2 = (p 2 )/2m, 
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where k& is the Boltzmann constant and m the atomic 
mass, is plotted against A' 5 for different detunings. The 
value of A' 5 is expressed in terms of the recoil energy 
gained by the atom after spontaneous emission, = 
h 2 k 2 x /2m. We note that for all detunings except A 5 = 
-40r, the results for \A' 5 \ > 125£r fall on a single line. 
We thus confirm the experimental finding in reference [15j 
that the depth of the potential generated by the 4 — > 
5 transition alone provides the appropriate scaling for 
large potential depths. The results for A§ = — 40-T have 
a slightly different slope. Since for this detuning the ex- 
perimental data in reference |15j only extends up to A' 5 ~ 
130£^r it is not possible to say whether this different slope 
is an experimental reality. 
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Fig. 1. One-dimensional kinetic temperature as a function of 
potential depth for different detunings (see legend). Simula- 
tions including both the F e = 4 and F e = 5 excited states. 



In figurc[2]wc compare to results of simulations includ- 
ing only the single excited state F c = 5. The most striking 
difference is that in figure [2] the results are essentially in- 
dependent of detuning over the whole range of potential 
depths considered. When the additional excited state is 
included there is a very clear dependence on detuning for 
shallow potentials, with the linear dependence extending 
further for larger detunings, giving rise to lower minimum 
temperatures. For larger detunings the linear behaviour 
is preserved to a lower potential depth, which makes it 
possible to reach a lower temperature. This phenomenon 
has been observed experimentally [T4 1 IT5 1 IT6] . and has up 
to now been at variance with all theoretical simulations, 
semiclassical or fully quantum mechanical. Thus, we can 
conclude that it is the additional excited state that causes 
this dependence of the point of decrochage on detuning. 

In figure [3] we show the minimum temperature T m i n 
achieved at different detunings, while figure 5] gives the 
laser irradiance at which this minimum was obtained for 
the same detunings. Around the minimum the simulations 
give a considerable amount statistical noise. This is partic- 
ularly true below decrochage, where a very small number 
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Fig. 2. One-dimensional kinetic temperature as a function of 
potential depth for different detunings (see legend). Simula- 
tions including only the single F c = 5 excited state. 



of atoms with very high momenta have a significant im- 
pact on the value of (p 2 )- In order to find the minimum we 
therefore made simulations for A' 5 in steps of 2Er around 
the minimum, and fitted the results to the functional form 



a\A' 5 \ + 6exp(- 



■c\A' 5 \ 



with a, b and c fit parameters. It 



was found that this function provides a good fit to simu- 
lated data to within the statistical uncertainties. We find 
that both T m ; n and the optimal potential depth follow a 
linear dependence on detuning. For T min this is consistent 
with the results in [14] over the range of detunings con- 
sidered. In [X4lJT6| it was also found that the minimum 
temperature is achieved for an optimal laser irradiance 
Zopt independent of detuning. Our results for I opt are dis- 
played in figure |4j The optimum irradiance varies almost 
a factor 2 over the range of detunings investigated, and 
thus clearly deviates from the experimental result. 




Fig. 3. Minimum temperature T m i n achieved for different de- 
tunings As, together with a linear fit. 
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Fig. 4. Laser irradiance 7 op t for which the minimum temper- 
ature is achieved as a function of detuning A5. 



4 Discussion 

At potential depths well above decrochage, the F c = 4 
level has very little influence. In reference [T5] it was rea- 
soned that for trapped atoms the dynamics is mainly de- 
termined by the lowest adiabatic potential, as most atoms 
get optically pumped into the extreme Mp sublevels. This 
potential is not affected by the F g — 4 — > F c = 4 transi- 
tion. If only the F g = 4 — > F c = 4 transition is considered 
the lowest adiabatic state has vanishing energy at all po- 
sitions (for detunings to the blue of the line) , and is thus 
dark to the laser light. This readily explains why, as long 
as the atoms are trapped in the lowest adiabatic poten- 
tial, their dynamics is determined only by the parameters 
of the F g = 4 — > F c = 5 transition. Even for large detun- 
ings the well known scaling of temperature proportional 



to I/A 5 oc A' 5 holds. 

Recently the dynamics of laser cooling in optical lat- 
tices was interpreted in terms of a bimodal momentum 
distribution [llII[2"tJll2"l"ll2"2] . The atoms are either in an un- 
trapped hot mode or in a cold mode where the atoms are 
trapped around a potential minimum. As atoms are trans- 
ferred from the hot to the cold mode the average kinetic 
temperature decreases. When the system is in steady state 
the interchange of atoms between the two modes is in bal- 
ance. For large potential depths essentially all atoms are 
trapped, giving rise to a truncated Gaussian momentum 
profile. At lower potential depths the hot mode can be ob- 
served even in steady state, giving rise to a deviation from 
a Gaussian velocity profile in the wings of the distribution. 

The dependence of the point of decrochage on detun- 
ing found in this paper and in reference [15] can also be 
understood from the bimodal picture of Sisyphus cooling. 
The cold mode is, as explained above, completely deter- 
mined by the F s = 4 — ► F e = 5 transition, and hence its 
temperature scales proportionally to A' 5 only. According 
to the bimodal model the hot mode starts to get populated 



0.2 0.3 
position z/X 

Fig. 5. The three lowest adiabatic potentials Ui(z) scaled by 
\A' 5 \ for As, = -Wr (upper panel) and As = -40T (lower 
panel). The solid line shows potentials calculated including 
both the F c — 4 and F c = 5 excited states, while the dashed 
line shows potentials calculated including the F c — 5 excited 
state only. In the upper panel the solid and dashed lines are 
almost identical. The potentials including only a single excited 
state scale with A' 5 only, and are hence identical in the upper 
and lower panel. 



around decrochage, thus driving up the value of {p 2 ), even 
though most atoms are still trapped in the cold mode [2 1 1 
122] . Even a relatively small population of the hot mode will 
dominate the value of {p 2 } since the atoms in this mode 
have no upper limit for their momenta, and for shallow 
potentials (p 2 ) may even diverge |23j . The increased pop- 
ulation of the hot mode is associated with atoms leaving 
the lowest adiabatic state. Our interpretation of the re- 
sults in figure [1] is that while the form of the cold mode 
is unaffected when the F g = 4 — > F c = 4 is taken into 
account, the rate of transfer of atoms from the cold to 
the hot mode, i.e. away from the lowest adiabatic state, 
is reduced. As the magnitude of this effect depends on 
the laser detuning from the F g = 4 — > F c = 4 transition 
a dependence of the point of decrochage on detuning is 
introduced in this way. 

In figure [5] we show the three lowest adiabatic po- 
tentials for detunings A5 = — 10r and A5 = — 40F. As 
noted above the lowest potential is identical for both de- 
tunings, giving the same dynamics for both cases. How- 
ever, as the potential depth is reduced the excited states 
in figure [5] gain significant populations. Since the poten- 
tials of these excited states are very different at different 
detunings the universal temperature dependence is vio- 
lated. The universal dependence persists to lower poten- 
tial depths at large detunings. We therefore conclude that 
the transfer of atoms to adiabatic states with higher en- 
ergies is more likely at small detunings, while for large 
detunings the potential depth has to be lowered even fur- 
ther before this transfer becomes important. As shown 
in figure E] at A5 = — 10 r (and for potentials where the 
F g = 4 — ► F e = 4 transition has been excluded) there are 
avoided crossings involving the lowest adiabatic potential 
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at z = A/8 and z = 3A/8 (A is the laser wavelength), 
while for A5 = —40-T there are distinct gaps. A tentative 
conclusion is therefore that the formation of this gap in- 
hibits transfer from the cold to the hot mode, although 
the details of this effect remains to be worked out. 
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Fig. 6. Momentum distributions for A$ = — 10.T (plus signs, 
black online) together with a fit to a double Gaussian, and 
for zis = —40-T (crosses, red online) together with a fit to a 
single Gaussian (double Gaussian would look the same). The 
momentum distribution has been binned into lpR-wide bins, 
where pr = Mr. The potential depth was IZI5 = 20Er, i.e., 
well below decrochage for A5 — —10r. 



This interpretation is also supported by the simulated 
momentum profiles. As an example momentum profiles 
at \A' 5 \ = 20E R at detunings A 5 = -lOT and A 5 = 
— A0r are displayed in figure El together with fits to dou- 
ble and single Gaussians respectively. (Fitting also the 
A 5 = — 40-T profile to a double Gaussian gives no im- 
provement as the widths of the two Gaussians in this case 
adjust to the same value, indicating that the distribu- 
tion really is well described by a single Gaussian.) For 
the smaller detuning the wider hot mode is clearly visi- 
ble. The fit gives for the cold mode (i.e., the central peak) 
widths corresponding to (p 2 ) = 1.9^ for A5 = — 40T 
and (p 2 ) = 1.4p| for A 5 = — 10r. Adding the hot mode 
gives a total (p 2 ) = 142p^ for A5 = —10r, even though 
the integral of the Gaussian reveal that both modes con- 
tain roughly the same number of atoms (49% hot, 51% 
cold). Even at the larger detuning a very small number of 
hot atoms increases (p 2 ) to 4.1 p^. Considering the more 
than one order of magnitude difference in the overall (p 2 ) 
between the two detunings, we find that the width of the 
cold mode is remarkably similar, showing that indeed even 
well below decrochage there is a significant population of 
the cold mode, with characteristics largely independent of 
the detuning. 



5 Conclusions 

In summary, we showed that at low potential depths, around 
the so-called point of decrochage, the temperature achieved 
by Sisyphus cooling does depend on both the potential 
depth and the detuning from resonance. For these po- 
tential depths it is necessary to include several excited 
hyperfine state in the theoretical description, in order to 
get accurate results. Simulations including only a single 
excited hyperfine state show no dependence of the tem- 
perature on detuning. This finding agrees very well with 
the experimental results of references [I4TITMT5] . previ- 
ously unreproduced by simulations. At larger potential 
depths, where the temperature depends linearly on poten- 
tial depth, we find that the additional excited hyperfine 
state has no effect. This is also in agreement with the ex- 
perimental results in reference |15| that the temperature 
scales with the potential determined including only the 
Fv = 4 — ► F e = 5 transition. 

We thank Anders Kastberg, Stefan Petra, and Peder Sjolund 
for many valuable discussions. This work was supported by the 
Swedish Research Council (VR) and by the EPSRC through 
grant number EP/D069785/1. 
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